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1 Introduction 

In this paper we present an explicit solution of the elliptic generalization of a quantum 
many body model in one dimension which is usually associated with the names of Calogero 
and Sutherland [Cll ISult IC2] and whose quantum integrability was proved by Olshanetsky 
and Perelomov |0P] . Our solution is based on a remarkable functional identity stated 
in Lemma 12.1.11 below. This identity was found in |Lll IL2] using quantum field theory 
techniques |CL] . but to make the present paper self-contained we also include an alternative, 
elementary proof of this key result. Our main results are explicit formulas for the eigenvalues 
and eigenfunctions of this model as infinite series. We also prove absolute convergence of 
these series in certain special cases. 



1.1 Background 

The elliptic Calogero-Sutherland (eCS) model is defined by the differential operator 

j=l j l<j<k<N 

^ langmann@kth . se 
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with —71 < Xj < 71 coordinates on the circle, A^ = 2, 3,..., 7 > —1/2, and the function 

V{r) = A ■ 2u \ .o w,i = Pir I vr, i/3/2) + co (/5 > 0) (1.2) 
4sin r + i jm) 2\ 

with the Weierstrass elhptic function p and the constant 

'° = 12 ~ 5i 2 sinh2[(/5m)/2] = 12 ~ ^ (1 - q^-f ^^'^^ 



where 



g = e-^/' (1.4) 

(see Appendix lA.ll for more details); we find it convenient to shift p by a constant to 
simplify some formulas later on. We also introduce the function 

00 

9{r) = sin(r/2) Y[{1 - 2?^" cos(r) + g*") (1.5) 

n=l 

which is equal, up to a multiplicative constant, to the Jacobi Theta function d\{r /2). This 
allows us to write ^ 

^(^) = -^log^W (1-6) 

(see Appendix lA . 1 1 for details). This differential operator H defines a quantum mechanical 
model of N identical particles moving on a circle of length 27r and interacting with a two 
body potential '-yVir) with 7 the coupling constant. To be more precise: the model we are 
interested in corresponds to a particularly 'nice' self-adjoint extension of this differential 
operator |KTj (which is essentially the Friedrichs extension [RSj ) . and we only consider 
eigenfunctions describing non-distinguishable particles. It is convenient to parametrize the 
coupling constant as follows, 

7 = 2A(A-1) (1.7) 

where A > 1/2 is the coupling parameter. Our restrictions on the parameters q and A 
are for simplicity and since they include most of the cases where H defines a self-adjoint 
operator bounded from below, but many of our results can be analytically continued to 
other complex parameter values. 

In the limiting case g = (/3 ^ 00) we have V{r) = (1/4) sin~^(r/2), and the differential 
operator H in (11. ip reduces to the one defining a celebrated model solved a long time 
ago by Sutherland jSulj who found explicit formulas for all eigenvalues and an algorithm 
to construct the corresponding eigenfunctions of H. The exact eigenfunction \l/n of the 
Sutherland model can be labeled by partitions n, i.e. integer vectors n = {ni,n2, . . . , n^r) G 
Nq such that fii > n2 > ■ ■ ■ > ni<^ > 0, and they are of the form 



with 



v[/„(x) = J-„(z)vl/o(x), z/='e^^^ (1.8) 
v^o(x)= n ^(^fc-^.)' (1-9) 

^<j<k<N 
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and the J-Jy'L) symmetric polynomials with simple corresponding eigenvalues given in (12.131) 
below. The are known as Jack polynomials and are also of interest in combinatorics; 
see jAfcPl lSt]. 



It is interesting to note that, to get all eigenf unctions of interest in physics, one should 
multiply the r.h.s. of (11.81) by a factor e"^'^^^^"' ^^^^ with an arbitrary P eN corresponding 
to an additional possible center-of-mass motion. One can account for this by extending the 
definition of the Jack polynomials to all integer vectors n such that ni > n2 > ■ ■ ■ > 
(with un possibly negative) as follows. 



for all partitions A and P G N; the ca,p are some non-zero constants depending on the 
normalization of the J^. This equation actually holds true also for negative integers P; see 
e.g. [St]. It is common to ignore this point and restrict the n to partitions. 

Sutherland's solution method is based on the fact that, for g = 0, the differential 
operator H in (II. ip has an exact eigenstate \l/o(x) of a form as in (II. 9p . This is no longer true 
for g > |Su2j , and thus Sutherland's approach cannot be generalized to the elliptic case. In 
this paper we elaborate another algorithm which allows to solve also the general elliptic case 
[Lit IL5j . In the trigonometric limit g = 0, this algorithm simplifies to one which is different 
from Sutherland's; see [L3] for a comparison of these algorithms. The crucial difference 
is that, while Sutherland computes the Jack polynomials jrn(z) as linear combinations of 

the following basis in the space of symmetric polynomials, M;v(z) = IljLi "^^^ 
partitions A and the sum over all distinct permutations tt in the permutation group S^^ we 
use a different generating set of functions /n(z); see (12.90 for 9(z) = {1 — z). It is important 
to note that the /„ are well-defined and non-zero also for integer vectors n which are not 
partitions. Thus the /„ provide an overcomplete generating set in the space of symmetric 
polynomials. Expanding the in this generating set 



the expansion coefficients a-a{™) obey simple recursion relations |L3j which can be solved 
explicitly |L5j . 

We emphasize that, in our approach, we get eigenf unctions \E'n for all integer vectors n. 
Obviously, if E-^ for a non-partition n G is different from all £\^pe for partitions A and 
P G N, then the corresponding eigenfunction \I'n must vanish, otherwise it must be a linear 
combination of the ^\-pe such that 8^ = ^^a-PgH It would be interesting to investigate 
this overcompleteness. 

1.2 The nature of our solution 

For g > 0, our algorithm leads to eigenfunctions \&n as in (I1.8p - (ll.l0p with a natural 
generalization of the functions /n(z) given in Proposition [2211] below. It is natural to regard 

^We have checked this for many different integer vectors n and parameters A in the case g = and 
TV = 2 using MAPLE. 



JA-Pe(z) = Cx,p{ZiZ2 ■ ' ' Zn) ^ J\{t) , 6 = (1, 1, ■ ■ ■ 



1) 




(1.10) 



m 
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the functions JJy'L) thus defined as an elhptic generahzation of the Jack polynomials. It 
is important to note that the latter are no longer polynomials for g > 0; see Remark 1 in 
Section [5] for a more precise characterization. Moreover, the corresponding eigenvalues E-^ 
and the coefficients an(in) are more complicated (LIJ. 

One main result in this paper is an implicit equation determining the eigenvalues and 
an explicit formula for the coefficients aj(m) depending on (Theorem I4.1.2p . We also 
present a solution of these equations by some variant of perturbation theory to all orders. 
This leads to fully explicit formulas for and aj(m) as infinite series (Theorem I4.3.ip . 
Our series are (essentially) power series in the parameter 7, but it is important to note that 
the small parameter is still and 7 is only a convenient book keeping device. A simple 
analogue illustrating this latter feature of our solution is the following function, 

with 5 the Kronecker delta: this function has a complicated power series in while its 
power series /(g) = X]^2 7 simple and can be written explicitly. Moreover, 

it is not difficult to prove that /, = 0(g2rfl) and that this series has a finite radius of 
convergence in . 

We also prove that our series for E-a, and aj(m) have a finite radius of convergence in g^ 
under a certain hypothesis which is fulfilled in several interesting cases, including N = 2 for 
non-integer A (Lame case), and the groundstate n = for all > 2 and irrational A. Our 
estimates to prove convergences are not optimal, and we believe that our series converge 
for all possible parameter values. In numerical computations it might be more efficient to 
solve the implicit equation for by iteration rather than using our explicit series. 



1.3 Previous results 

We now discuss previous results related to our work. 

For N = 2 the eigenvalue equation of the eCS differential operator can be reduced to 
the so-called Lame equation which was studied extensively at the end of the 19*'' century 
(see |WWj for a summary of these classical results) and more recently in {EK\ IR2j . It is 
known that, for integer values of A, the Lame equation has a finite number (depending on 
A) of eigenf unctions and corresponding eigenvalues which can be computed algebraically 
|WW] , and a generalization of this result to > 2 was found by Gomez-UUate et. al |GGR] 
and led to particular many-body generalizations of the Lame equation which, however, is 
different from the eCS model. There exist various other interesting results on the solution 
of the eCS model for integer values of A: Dittrich and Inozemtsev obtained explicit formulas 
for the eigenvalues of the eCS model for the cases A = 2 and A^ = 3 |DIt HI] and A = 2 and 
general A^ [12]. A Bethe ansatz solution of the integer-A eCS models and, more generally, 
the elliptic Ruijsenaars model |Rlj (which reduces to the eCS model in a certain limit), 
was given by Felder and Varchenko |FV2l IFV3j : see also [HI IHSYt [T] . There exists an 
interesting separation-of-variable approach to the solution of the A^ = 3 eCS model by 
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Sklyanin [Sj. The latter approach seems to be similar to ours in that it is also based on an 
identity closely related to the one in Lemma 12.1.11 and it was elaborated in detail in the 
trigonometric limit g = by Kuznetsov et.al. |KMS] . We also mention the work of Niifiez 
et.al |FGP] who computed the eigenvalues of the eCS model for arbitrary values of A up 
to 0{q'^) for < 4 and up to 0{q'^) ioi N = 2 and which our results extend to arbitrary 
orders and particle numbers. Other related work is by Komori and Takemura |KT] who 
presented a perturbative algorithm to solve the eCS model and proved that Schrodinger 
perturbation theory has a finite radius of convergence in g^; see also As far as we can 
see is the approach in [KT] different from ours (we obtain results which are different from 
standard ones even for g = |L3j ) and has not yielded results equally explicit as ours. A 
complimentary discussion of known results about eCS type systems and further references 
can be found in jKRj . 

As mentioned, our approach is based on a remarkable functional identity in Lemma P2.1.1[ 
This identity can be regarded as a natural quantum analogue of the Backlund transforma- 
tion discovered by Wojciechowski [W]; see also [KS]. Similar identities were obtained by 
Felder and Varchenko [FVlJ and, for = 3, by Sklyanin [S]. Moreover, its trigonometric 
limit is equivalent to a well-known identity for the Jack polynomials due to Stanley; see 
Proposition 2.1 in Ref. ^tj. Similar identities exist for a large class of Calogero-Sutherland 
type models and provide a tool to construct explicit formulas for the eigenfunctions of all 
these systems in a unified way [HLj . It is interesting to note that the functions /n(z) for 
g = are the building blocks of this explicit solution for all these models. 

The results presented in this paper evolved over several years. The key to our solution 
was found already in 2000 [Llj . but we realized only in 2004 that it is possible to obtain 
an explicit solution to all orders |L5j . Readers interested in how this result emerged can 
consult the two earlier versions of the present paper |L4j . The first version ([L4J vl) is 
a manuscript from 2001 which remained unpublished until January 2004. It contains our 
solution algorithm together with a recipe how to compute the eigenvalues as power series 
in the nome g of the elliptic functions. In the second version ( |L4j v2) published in April 
2004 we added explicit formulas for the eigenvalues up to O(g^) for = 2 and up to O(g^) 
for < 4. These formulas were obtained by straightforward but tedious computations. 
These latter results suggested an alternative, more efficient, method allowing to find the 
solution to all orders in g^. This result was previously announced in Refs. |L5l IL7j and is 
elaborated and extended in the present paper. 

A key issue in this evolution of our result was that we realized that there exist better 
answers to the question: "What parameter should we use to expand our solution in?", 
than the obvious one. The obvious answer is: "g" |Llj . However, this leads to very 
complicated formulas already at low orders. A better answer is: "7", since it leads to 
much simpler formulas at low order; see Section [3l The final answer which allowed us 
to deduce explicit formulas to all orders is: "77", where t] an auxiliary parameter which 
we introduce into the problem and which can be set to 1 at the end of the computation. 
In particular, the explicit formulas for the eigenvalues which we thus obtain are of the 
form = ^o(n) + Yl^=i'U^^rn{p) where £mip) = C(g^™) are Taylor series in g^; see 
Theorem 14.3.11 and Remark I4.3.2[ 
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1.4 Plan of the rest of the paper 

In the next section we summarize preliminary results which we need. In Section E] we derive 
the solution of the eCS model, i.e. its eigenfunctions and corresponding eigenvalues, as a 
series in 7 up to order 0(7^) and 0(7^), respectively, using an elementary argument. Our 
main results are obtained in Section HI explicit formulas for the solution of the eCS model 
as a series in 7 to all orders, together with a sufficient condition for absolute convergence 
( Theorems I4.1.2l and l4.3.1]) . We end with remarks in Section O Details of our computations 
and some proofs are deferred to three appendices. 

Notation: We denote as Z, Z', Nq and N the sets of all, all non-zero, all non-negative, and 
all positive integers, and C and M are the complex and real numbers, respectively. We use 
bold symbols for vectors with components, e.g. n G is short for {ni, . . . ,71^) with 
rij G Z etc. The symbol "5" always means the Kronecker delta, in particular, 6{m, n) = 
rijLi ^rnj,nj for m, n G Z^. For x G M we denote as \x~\ the smallest integer larger or equal 
to X. Definitions are indicated by the symbol "=" . 



2 Summary of preliminary results 

In this Section we collect a few results which we need. 
2.1 A remarkable identity 

The starting point for our solution method is the following. 
Lemma 2.1.1 Let 

J-,/ ^ def Y[l<j<k<N^i^k — Xj) Y[l<j<k<N 

with 6{r) in U.5\) . where x = {xi, . . . ,xn) G and similarly for y . Then the following 
identity holds true, 

E(^-9T)^(^'y) = 2MA-l) {v{xk-x,)-Viy,-yk))Fi^;y) (2.2) 

j=l l<j<k<N 

with V{r) as in ( fi.6j) . 

(Proof in Appendix IA.3I ) 

As already mentioned, the result in Lemma 12.1.11 was obtain in |L2j using quantum 
field theory techniques. The proof given in Appendix IA.3I is elementary and based on the 
following functional identity [WWj 

(j){x)(j){y) + (j){x)(j){z) + (j){y)(j){z) = f{x) + f{y) + f{z) ii x + y + z = (2.3) 
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where 

0(x) = ^\oge{x) , fix) = '^[V{x) - - Co] (2.4) 
and Co the constant in (11. 3p : a proof of the latter identity is in Appendix IA.2[ 
It is important to note that one can write the identity in (12.21) as follows, 

i/(x)F(x;y) = /7(y)F(x;y) (2.5) 

where H is the differential operator in (II. ip but acting on different arguments x and y, as 
indicated. 



2.2 Reformulation of the eigenvalue problem 

From Lemma 12.1.11 a straightforward computation leads to a result which is the next step 
in our solution. To state this result we find it convenient to define 

{E,k)e = Sje-6M Vj, fc, £ = 1, . . . , iV, j<k (2.6) 

and use the shorthand notation 

£> l<j<k<N u&Z 

for functions a on integer vectors in Z^. 
Proposition 2.2.1 Let 

^n(x) =7n(z)^o(x), u & and Zj = e'""^ (2.8) 
with \E'o(x) defined in ^1.9\) and U.5\) and the special functions 

with the integration paths 

Cj : = e^'e''^' , -tt < (/^j < vr , < ei < 62 < ■ ■ ■ < En < (3 (2.10) 

and 

00 

w = (1 - ^) n [(1 - ^""'') (1 - ^""/^)] • (2.11) 

m=l 

Then the eCS differential operator H in U.l\) - in~^) obeys 

iJF„(x) = ^o(n)F„(x) - 7 5^ 5,F„+£,(x) (2.12) 
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where 

N 

So{n) J2 (^^- + + 1 - 2j)Y (2.13) 

and 

S. = W\YZ^ Vz/gZ', ^0 = 0. (2.14) 

(Proof in Appendix iBl) 

This result concluded our discussion in [L2j Jl Our proof in Appendix [B] is (essentially) 
by expanding the identity in (12.51) in a Laurent series in the variables C,j = e'^J in a certain 
region in and equating the expansion coefficients. 

As discussed below, the integrals in fl2.9p -( l2rTT]) are well-defined and independent of the 
parameters ej in the specified range due to Cauchy's theorem. It is important to note that 
Su = for z/ < and g = 0. As we will see, this simplifies the solution for g = drastically. 

An important consequence of Proposition 12.2.11 is the following. 
Corollary 2.2.2 Let 

^(x) = J2 «(m)F^(x) (2.15) 

with coefficients a{m) satisfying the following relations, 

[£o{ui) - £]a{m) = 7(§a)(m) (2.16) 

for some constant S and 

(§a)(m) J2 - »>) (2.17) 

with Fn(x), SqIii), and S^, defined in Proposition \2.2.1\ and u in ^2. 7\ ). Then \E'(x) is an 
eigenfunction of the eCS Hamiltonian in U.l\) with eigenvalue S: /7\E'(x) = £^\E'(x). 

This is a simple consequence of the ansatz (12.151) and Proposition 12.2.11 which imply 
{H - = 5^([£:o(m) - £]a{in) - 7(§a)(m))F^(x) . 

m 

Thus the problem of solving the eCS model is reduced to finding solutions S and a(m) of 
2.3 Properties of the functions /n(z) 

The functions /n(z) defined in (I2.9p - (l2.1ip play an important role in our solution. To see 
that they are well-defined we note that 

< e{\z\) < \e{z)\ < e{-\z\) < oo for <\z\<l (2.18) 
■^Note that /n(z) here was denoted as ■p(n;x) in [L2| . 
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(see Appendix IC.ip . and thus the integrand in (12. 9p is analytic in the region 1 < |^i| < 
1^2! < ■■■ < \C,n\ < if all \zj\ = 1. Thus the integral in (12. 9p are well-defined and 
independent of the closed integration paths Cj as long as they are within this region of 
analyticity (Cauchy's theorem), which is obviously the case for the ones in (l2.1Up . 

It is important to note that the /n(z) are symmetric functions: fn{zi, Z2, ■ ■ ■ , Zn) = 
/n(-27r(i), -27r(2), • • • , ^7r(7V)) for all 71 G Sn- Morcovcr, they are polynomials for g = |L3j . 
As already mentioned, for q = the /n(z) provide a basis in the space of symmetric 
polynomials if the n are restricted to partitions, even though they are, in general, non-zero 
even for non-partition integer vectors n \L3\ IHL] . We expect that similar results hold true 
also for g 7^ 0. 

We finally state an upper bound for the functions /n(z) which we need to prove square 
integrability of our eigenfunctions: 

Lemma 2.3.1 The functions /n(z) defined in Ii2. 9\) obey 

|/„(z)| < Cg^^(^l"^l-^^"^") , \zj\ = 1 (2.19) 

with the constants 

j^mJL_ k'^^ c '-^ [2e(-g^)]^(^-^)v^ 

N + h' 1 + 26' [(1 - g^~^)e(g2-2fe^)/(l - g2-26K^]7V2A ^ • ^ 

for arbitrary b > 0. 

(Proof in Appendix IC.li ) 

Note that < K < K/2, 2 — 2bK > 0, and C < oo. It is worth mentioning that these 
estimates are not optimal (and in fact meaningless for g = 0), but they are good enough 
for our purposes. 



2.4 Lagrange's reversion theorem 

We finally recall a well-known result due to Lagrange which will play a key role in our 
solution: 

Theorem 2.4.1 (A) Let (f{z) and g{z) be functions of the complex variable z analytic on 
and inside a closed contour C surrounding a point z = a, and rj a complex number such that 

\rnp(z) — a\ , , 

sup ' ^ , ^ ^ , ' < 1 . (2.21) 

zee \z-a\ 

Then the equation z = r]{p{z) has a unique solution z = C, inside of C, g{C,) can be expanded 



as 



m=l 

and this series is absolutely convergent. 



,(0 =«W + I:;L_ (,(„).._,(„) j. (2.22) 
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(B) The result in Ii2. 22\] holds true in the sense of formal power series even for values of rj 
where Ii2.21\) is not fulfilled. 



Proof: Part (A) is equivalent to Lagrange's theorem as stated in [WW], Section 7.32. The 
result in (B) is a simple implication. □ 

Remark 2.4.2 We note that Ii2.22\) does not rely on analyticity and applies to even more 
general situations where one knows that (p{z) and g{z) have formal power series in [z — a) 
hut has no information about convergence of these series. Indeed, inserting 

oo 

if{z)=Y,Mz-<^T (2.23) 

n=0 

into the equation z = rnp{z) and making an ansatz ^ = a + 'Y2'^=i^kV^ /^^ ^^■^ solution one 
can solve for the recursively and compute 

^ = a + ?7V9o + r^VoV^i + V^i^l^2 + Vovl) + v'^i.vl'f?. + 'ivlviV2 + Vov\) + ... ■ (2.24) 
Inserting this in g{z) = X]-^o5'"('^ ~ '^)" yields 

+ 9iyi'^2 + f^o^?]) + t1^{9avI + 3^3^o<^i + ^2[2<^o<^2 + 3(^o<^i] 

+ S'lboV^s + 3v3oV5iV52 + V^oV^?]) + • • • (2.25) 

which can he straightforwardly extended to higher powers in rj and agrees with Ili2.22\) to all 
orders in rj. 



3 Perturbative solution 

We will present our solution to all orders in the next section. As a heuristic motivation, 
we present in this section a pedestrian approach allowing to compute the first few terms of 
our series solution in a simple manner. 

We found that (I2.16p can be efficiently solved by making the ansatz 

oo oo 

a(m) = J]7^a(^)(m), £ = ^^'£^'^ (3.1) 

leading to the following system of equations equivalent to fl2.16p . 

[^o(m)-^(°)]a(°)(m) = (3.2) 

and 

s 

[£o(m) -^(°)]a(^nm) = ^£'^''^a'^'~''\m) + {^a^'-^^){m) Vs > . (3.3) 

s'=l 
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As mentioned in the introduction, 7 here is not assumed to be small but only serves as a 
convenient book keeping parameter to organize our solution. It is important to note that, 
for each n G Z^, the initial condition in (13.21) has the following solution, 

4°)=^o(n), «(°)(m) = 5„(m) (3.4) 

where (5n(m) = 6{m, n). Note that an^°^(n) is then undetermined by (13.31) and could be set 
to any value. This ambiguity corresponds to the freedom of multiplying the eigenfunction 
by an 7-dependent constant. Our choice below is a convenient normalization. 

Equation (13.31) implies 

£i'^ = -{Sa^^-^^)in) Vs > (3.5) 
where we set an''(n) = for all s > 0, and 

a('^)(m) = -— i -JY.^n^^t''\^) + (S«(^-^))(m)) V. > ; (3.6) 

6n(m-n)V^ / 



we use the convenient shorthand notation 



if m = n 



bn{m — n) ]^ [i^o(m) — £^o(n)] ^ otherwise 



(3.7) 



We thus can compute, recursively, all Sn^ and an\m). We obtain 



bJr>) 



(5(m - n, i>i + U2) 



al'^m) 



V 9 9 - 



P(3) _ _ \^ Q C C ^(0,t>l + t>2 + t>3) 



«„lmj- 2^ ^.,^.,^.3 1^ ^(i>i)6n(i>3)2 



5(m - n, i>i + i>2 + i>3) 



C(4) _ C C C C f 

t-'n / ^ >~'Pi>~>V2^Vi>~>V4,\ 



5(0, l>i + i>2)5(0, I>3 + i>4) 



^n(i>l)&n(l>3)2 
5(0, l>i + l>2 + l>3 + Ua) 



+ i>2 + ^'3)&n('>l + i>2)bnil>l] 



(3. 
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etc. It is straightforward but tedious to continue this computation to higher orders in 7. 
Below we present a more efficient computation method allowing us to derive closed formulas 
for £^n and an(m) to all orders in 7. 

It is straightforward to expand ^n'''' and in powers of g^, but the resulting formulas 
are rather complicated, and we therefore do not present them here0 

We thus find, for each n G Z^, an eigenfunction ^'n(x) and a corresponding eigenvalue 
as a formal power series in 7. However, it is important to note the following potential 
problem: It is possible that So{in) — So{n) vanishes for some m 7^ n with Ylji^j ~ 
which appear in the above sums (recall that these series are built of terms 6n(m ~ n) in 
(13. 7p with m = n + X]k=i ^^"^ such m obey this latter condition), and in this 

case the sums are not defined. We refer to such m as resonances. Thus the results in the 
present section make sense only if there are no resonances, i.e. if 



Since there exist parameters A and n where resonances exist this is an important restriction. 

We now discuss cases with resonances. For example, for = 2 we have 

£^o(ni) - £^o(n) = 2z/(z/ + ni - 71,2 + A) , z/ = mi - = -(m2 - ^2) (A^ = 2) , (3.10) 

and thus there is a resonance for v = n2—ni — A provided that A is an integer. It is easy to see 
that there exist resonances for all n and iV if A is an integer, but for N > 2 there exist even 
resonances for non- integer A. For example, for = 3 and n such that u = (rii — 2?t,2 + n3)/3 
is a non-zero integer, there is a resonance at m = (^i — z/, n2 + 2z/, n-^ — u) for arbitrary values 
of A. Such A-independent resonances exist for all > 2. It is also interesting to note that, 
if A^ > 2, there exist infinitely many ways to represent a resonance as m = n + Yl^=i 
for m > 2: this is true since 



^hh + E,2,3 + • • • + E,,_,,, - E,,,, =0 if 1 < Ji < J2 < • ■ ■ < < AT, / < AT . (3.11) 



Due to this identity resonance can occur in the above sums infinitely many times. 

However, it is important to note that there exist cases without resonances where the 
results of the present section apply. In particular, (I3.10p shows that no resonance exists 
for N = 2 and non-integer A. Moreover, results from numerical experiments performed 
with MAPLE suggest that, for arbitrary particle numbers and non- integer A, there exist 
infinitely many n without resonances. 

Anyway, we show in the next section that this potential problem of resonances can be 
circumvented by resummations. 

4 Solution to all orders 

We now introduce some notation which allows us to extend the solution above to all orders. 
''The interested reader can find a few such formulas |L4j v2. 



N 




(3.9) 
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4.1 Implicit solution 



We write fl2T^ as 

(A-£)a = 7§Q; (4.1) 

with 

(Aa)(m) =fo(m)a(m) . (4.2) 

It is natural to interpret A and § as linear operators on the vector space Map(Z^;C) 
of functions a : ^ C, m i— a(m). (The topology of this space plays no role in 
our discussion and is therefore ignored.) It is also convenient to introduce the following 
projection on Map(Z^; C), 

(P„a) (m) = 5(m, n)a(m) . (4.3) 
The equation in (14. ip can be solved using the following simple but powerful result. 

Lemma 4.1.1 Let A and B he linear operators on a vector space V and consider the 
eigenvalue equation 

{A-£)a = Ma (4.4) 

with the eigenvalue £. Let P he some projection on V commuting with A and = / — P. 
Then 

a = [J + (A-£ -P^B)-ip^B]ao (4.5) 

is a solution of this eigenvalue equation provided that £ and satisfy the following condi- 
tions, 

Pao = ao (4.6) 

and 

£a^ = Aao - PB[/ + (A - ^ - P^B)-¥^B]ao . (4.7) 
Proof: Applying the projections P to (14. 4p and inserting 

a = ao + P"'"a, ao = Pa 

we obtain 

£ao = Aao - ™a (4.8) 
where we used that A and P commute. In a similar manner, applying P-*" to (14.41) . we get 

AP^a - £F^a - P^BP^a = P^ 



Obviously the last two equations together are equivalent to (14.41) . From the last equation 
we get 

F^a = {A-£- P^B)-¥^Bao , 
which implies (14. 5p . Inserting (14. 5 p in (14.80 we obtain the condition in (14.70 . □ 

If we apply this result to (14.11) using the projection P = Pn and B = 7§ we can compute 
a solution a = a^ and £ = £n- It is easy to solve (14.61) : its general solution is ao = c5n 
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with 5n(m) = 5(n, m) and an arbitrary constant c. It is obvious that c is an overall 
normalization constant which can be set to 1 without loss of generality. This is equivalent 
to our normalization condition an\'n.) = 1 and = for s > 1 in the previous section. 
We thus obtain 

ttn = [/ + 7(A - ^ - 7P^§)-^P^§]5n , 

and the condition in (14.81) determining the eigenvalue becomes 

£n = Soin) - 7(§«n)(n) . 
We expand in a geometric series 



oo 

s=0 



where we introduce the convenient notation 

1 



[[^o(m)-^]]„ 
which allows us to write 

{[A- S]-^F^a){m 

We thus obtain our first main result. 
Theorem 4.1.2 Let n e and 

^n(x) = 

where 



if m = n 

[So{m) — S]^'^ otherwise 



-a(m) . 



[[^o(m) - S]l 



a„(m)Frr,(x.) 



a 



s=i />i,...,i>, K=i nK=i ^^o(n + xir=i - 



(4.9) 



(4.10) 



(4.11) 



with -Fn(x), Soip), Sy, i> defined in Proposition \2.2.1[ and Sn a solution of the following 
equation, 



s=2 



(4.12) 



Then \l'n(x) is an eigenfunction of the eCS differential operator H in with 
eigenvalue S^: if\l/n(x) = £^n^n(x). 
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Our arguments above prove that this result holds true at least in the sense of formal 
power series in the nome q of the elliptic functions, but, as we discuss in the next section, 
in many cases the series have a finite radius of convergence. 



As before can the sum in (14.101) be restricted to m G Z satisfying Ylfj=i{^j ~ ^j) = 0- 



Note that (14.121) is an implicit equation determining the eigenvalues of the eCS model. It 
is straightforward to obtain from it a solution as a formal power series in 7 and thus recover 
the results in the previous section. It is important to note that the problem of resonances 
has disappeared in (14.111) and (I4.12p . and it only reemerges if one solves these equations by 
a particular series expansion (which is essentially the one discussed in the previous section) . 
However, it is possible to compute other series solutions avoiding resonances. 



Remark 4.1.3 It is important to note that Iji4-1^ simplifies in the trigonometric limit as 
follows, 

£n = £oin) forq = (4.13) 

(since for q = Q, Sy is non-zero only for z/ > 0, hut all the Kronecker deltas in the sums on 
the r.h.s. of Iji4-i^ (^'^^ ^ero if all Ur > 0), and we thus recover the well-known eigenvalues 
of the Sutherland model ISulf . Thus in this case, Theorem \4.1. ^ provides fully explicit and 
well-defined eigenfunctions \E'n(x) one can prove that for all partitions n, ^^o(n+X]^=o " 
£^o(n) > for all k > and > j.L3j . and thus resonances do not appear in \4.11\) for 
q = 0. As already mentioned, for partitions n and for q = these eigenfunctions are the 
same as Sutherland's IL3[ \HL^ . In the following we assume q > 0. 



We now introduce some useful notation. We write (14.121) as follows, 

£n = $n(^n) where Sn = £n " ^^o(n) 



and 



s=2 



(4.14) 



(4.15) 



'^=1 n^=i [[^0 (n + ELi i>i)-£oH- 
is a complex valued function of one complex variable z. Similarly we write (14. lip as 

a„(m) = G'n(4; m) (4.16) 

where 

5(m,n + ^^^^ j>«) 



G„(^;m)^='-5(n,m) + ^7^ J2 H^-: 



s=i ,>!,...,£<, K=i n«=i ^o(n + Er=i*^^) -^o(n) 
(for fixed n and m) is also a complex function. 



(4.17) 
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4.2 Analyticity results 



In this subsection we establish results concerning the analyticity of the functions $n(^) and 
Gn{z; m) introduced at the end of the previous section. This is rather technical but needed 
to prove convergence of our series solution. 

As we show below, for sufficiently small values of q, the above mentioned functions are 
analytic in certain regions of the complex z-plane. To simplify our discussion we now state 
a certain non-degeneracy condition which, when fulfilled, allows us to give a simple proof of 
convergence. We emphasis that neither this hypothesis nor the estimates we deduce from 
it are optimal, which is why we only get a proof of convergence in certain special cases. 

Hypothesis 4.2.1 The model parameters N, A and n G are such that there exist 
constants a G M and A > |a| satisfying 

N 

\£oim) - £^o(n) - a| > A > Vm ^ n such that "^{mj - Uj) = . (4.18) 
Note that a and A can depend on n, and A. 

It is interesting to note that, for rational values of A, it is trivial to find constants a ^ 
and A > |a| such that the condition in (14.181) is fulfilled: 

Lemma 4.2.2 Let A be rational and m and n > co-prime integers such that A = n/m. 
Then any of the following parameters 

a = ki + Xk2 + Oq, fci, A;2 € Z and A = |ao| < — (4.19) 



2m 



obey the condition in l[4.19^ . 



Proof: Equation (I2.13P implies that £^o(n) — is always of the form Vi + Az/2 for some 

integers Vi^2i and thus it is obvious that \vi + Az/2 — a| > A for all integers vi^2 and the 
parameters in (14.191) . □ 

This lemma is enough to establish that the functions defined at the end of Section HTTl are 
analytic in non-trivial z-regions for sufficiently small values of q (Proposition 14.2.31) . How- 
ever, to prove that the prerequisites for Lagrange's theorem are fulfilled (Corollary 14. 2. 41) we 
need that A is strictly larger than |a|, and it is this which is non-trivial in our hypothesis. 
At the end of this section we prove that Hypothesis 14.2.11 is fulfilled for all n if A^ = 2 
and A 7^ N (Lemma 14.2.61) . For A^ > 2 we only have some partial results, including the 
groundstate n = and irrational A (Lemma I4.2.7p . 

We now are ready to state our analyticity result: 

Proposition 4.2.3 Let a G M and A > be such that the condition in ( [^.igp holds true, 
b > a parameter which can be chosen arbitrarily, and q sufficiently small so that 



^ - g2/(7V+fe))3 (4.20) 
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satisfies B < A. Then ^n{z) in ( [^.i5p and Gn{z;in) in i4.11\) are analytic functions in 
the following region 

\z-a\</\-B (4.21) 
of the complex z-plane. Moreover, in that region the following estimates hold true, 

< A-/-|.-a| (^-^2' 

and 

\Gn{z;in)\ < 5(m,n) + (4.23) 

A — B — \z — a\ 

with the constant K in ^2.20) . 

Proof: The estimate in (14.181) implies 

1 ^ 1 ^1 

£o(jn) — Eoin) — z ~ \£q(u)) — £q(jli) — a\ — \z — a\ ~ A — 1^; — a| 

for all m = n + 7^ n and I2; — a| < A, and thus (I4.15P yields 



00 , X s-l 



=2 

where 

s 

Ks{^) = \lY E n^-^(^'^-i^«)>0 (4-24) 

£'i,...,£'s K=l 

for all m G and s G N; we used S*!^ > 0. One can prove that 

if,(m) < g2E, j-./(JV+&)5« (4.25) 

with 5 in (I4.20p and 6 > arbitrary; see Appendix IC.2I for details. This implies 

s-l 



provided that \z — a\ < A — B. Summing up the geometric series we obtain the estimate 
in (I4.22p . It is obvious that each term in the series defining $n(^) in (I4.15P is an analytical 
function in the region 1^ — a| < A, and thus our estimates above prove that ^n{z) converges 
and defines an analytic function in the region defined in (I4.2ip : see e.g. [WWj . Section 5.3. 

The proof of analyticity of G'n(-2; m) in (14.170 is similar: we now can estimate 
\G^{z- m)| < ^(n, m) + E ( _ j Ks{m-n) 

s=l ^ ' ' ' 

with Kgim) in (I4.24p . Using the estimate in (I4.25P we obtain 

|G„(^;m)| <5(n,m) + yg2E,iK-%)/(A^+^') ( ^ Y 

^-^ \ A — \z — a\ I 
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which imphes the estimate in f l4.23p provided that \z — a\ < A — B . Similarly as above we 
conclude that Gn{z; m) is analytic in the region defined in fl4.2ip . □ 

We solve the equation in (14.141) by using Lagrange's theorem 12. 4.11 The following result 
gives a sufficient condition that the prerequisites for the stronger version (A) of this theorem 
are fulfilled: 

Corollary 4.2.4 Let a G M and A > \a\ be such that the condition in lji4-iS\ ) holds true, 
b > an arbitrary parameter, and q sufficiently small so that B in ^■20{ ) satisfies B < 
(A — |a|)/3. Then for any of the following closed contours in the complex z-plane, 

C: \z-a\ = ^(^A- B + \a\- e>y{A -B- \a\f - 45^) , -1 < £ < 1 (4.26) 
the following condition holds true, 

sup ' ^ , ' < 1 , 4.27 

\z - a\ 

and $n(-2) o-nd all functions Gn{z\m) are analytic inside and on C. 

Proof: It follows from the estimate in (14.221) that the condition in (I4.27P is implied by 

B^ 



A- B -\z-a\ 
The latter holds true if and only if 



+ a < \z — a\ . 



{\z-a\-{A-B + \a\)f < {A - B - \a\f - AB^ , 

which has non-trivial solutions as in (I4.26P provided that 2B < A — B — \a\. This proves 
(I4.27P under the given assumptions. The statements concerning analyticity are implied by 
Proposition KTE\ □ 



Remark 4.2.5 The best possible bound in is obviously obtained for b = 0, but we 



need b > to prove square integrability of the eigenf unctions. 

To show that Hypothesis 14.2.11 is relevant we now discuss a few special cases where it 
holds true. 

Lemma 4.2.6 Let N = 2 and A be non-integer. Then Hypothesis \4. 2. 1\ holds true with 

a = 0, A = min|2i/(i^ + ni-n2 + A)| >0 ifX^Z {N = 2) . (4.28) 

Proof: This is a simple consequence of Equation (I3.10p . □ 

It is interesting to note that A in (I4.28P can be arbitrarily large, e.g. for ni— n2+A = k+r, 
k eN and 0<r<l/2, A can be as large as 2rk. It thus is possible to find cases where A 
is large enough for our series to converge for all values g < 1. 

The reason why we did not find a simple, general proof of convergence are the resonances 
discussed at the end of Section [31 as can be seen by the following result: 
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Lemma 4.2.7 Let N >2 and X and n G Z such that the no-resonance condition in liS. 9\) 



is satisfied. Then Hypothesis \4-2. 1\ holds true for a = and some A > 0. 

Proof: Since the set of all m such that |£^o(m) — ^o(n)| < 1 (say) is obviously finite, 

A^='inf|^o(m)-^o(n)| 

m 

can be computed as minimum over a finite set, and thus fl3.9p implies A > 0. □ 

Note that the no-resonance condition is satisfied in all cases n where there is no A- 
independent resonance and where A is irrational. Another interesting special case where 
Lemma [4.2.71 applies is the groundstate n = and irrational A. 



Remark 4.2.8 The careful reader might wonder why our estimates in Proposition \4.2.3 
are not analytic in . The reason is that these estimates are not optimal. Indeed, since 



^n{,z) = for q = (see Remark \4.1.3\ ), we expect that ^n{z) should vanish like q^ as 



g — > 0. In Remark ] C. 2. 1\ (Appendix) we motivate the 

Conjecture: Ks(0) < B'q^^^^ with B ^(^ -J-)H _ .4 29) 

(1 — q'^i"Y 

This would imply the following improved estimate 
Conjecture: \^Jz)\ < Bq^ ( (^- k-^l)^ V- ^N-nr^ _ u _ ^1)" - 1 

analytic in q^ and consistent with $11(2) = 0{q^). We expect there exists a similar improved 
estimate for Gn{z\ m). 



4.3 Explicit solution 

From the implicit equation for in Theorem 14.1.21 it is straightforward to get an explicit 
expression by expanding in 7 and thus extend the series solution in Section [3] to all orders. 
It is important to note that this amounts to an expansion around £ = £^o(n). As will be 
shown below, it is possible to expand around any point S = So{n) + a, a arbitrary, and thus 
the 'resonance denominators' So{m) — £o{n) can be moved to So{m) — So{n) — a. In this 
way the resonance problem can be circumvented, and we obtain the following fully explicit 
result. 

Theorem 4.3.1 (A) Let n G and a G M such that the condition in ( [/^.igp holds true 
for some A > 0. Then the eigenvalue equation if\t'(x) = S^^^x.) of the eCS differential 
operator in U.l\) - n~^) has a solution S = Sn as follows, 

00 00 

£n = Soi^i) + a + ^ ^ SiYlT=o^r, m) 

m=i eo,ei,-,em-i=o 

r^('')/ \'\£ 

xSiEZirir, m - l)(m - 1)! n ^ /; ^ ^ (4-30) 
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with 



s=2 £'i,...,i>s K=l fci,fc2,.--i'Ss— 1=0 

5(0, EU^.) 



Ul=\ ^o(n + Er=i - ^o(n) 



(4.31) 



and the notation defined in ^2. 0^ . I[2. ?| j and l[4-9^ - Moreover, the coefficients an(m) giving 
the corresponding eigenfunction \E'(x) = \I/n(x) according to Theorem 4-i-2 are 

oo m—l oo 



a 



„(m) = GW(a;m) + E '^(ELo^o ^)5(E^=ir^n ^) 

m=l ^=0 eo,ii,...,ii=0 

X (m - l)!(m - £)G(r"^)(a; m) J] ^"^"^y^'" 



(4.32) 



r=0 



wzt/i $nHo) as above and 

oo s oo 

G^;\a-ni) = 6r,oS{ui,n) + Y,Y E H E '^ELi^^, 

>s r=l fci,fc2,.--,fc3=0 

(5(m,n + ELi 



s=l £'i,...,j>s r=l fci,fc2,.--,fc3=0 



Uk=i ^o(n + Er=i^^) -^o(x) 



(4.33) 



The results above hold true for arbitrary parameter values, and Sn and an(m) above are 
independent of the parameter a in the sense of formal power series. 

(B) If A > \a\ and if q is small enough that B in Iji4-^0\ l satisfies the condition B < 
(A — |a|)/3, then all series above converge absolutely, the following estimates holds true. 



\Sn - So{n) -a\<^(^A-B + \a\- ^ {A - B - \a\y - 4B^^ 



and 



|an(m)| < 5(m,n) + g^.^'^K-"^) 



2B 



A-B-\a\ + y/{A- B- \a\y-AB^ 



(4.34) 



(4.35) 



with K in ^2. 2(A) . and the eigenfunction ^E'n(x) determined by the formulas above is square 
integrable. 



Proof: We observe that (14.141) is of the form z = f]y^{z) if we identify 

z = 4 , VVi^) = '^n{z) . (4.36) 

We thus can use Lagrange's theorem 12.4. II to solve (14.141) . Note that t] can be regarded as 
a book keeping parameter which can be set to any value. We set = 1 in the following. 
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We first prove the result under the assumptions stated in (B). According to Corol- 
lary |4]2]4] the prerequisites for Lagrange's theorem are then fulfilled for any of the contours 
defined in ( 14.261) . Thus the equation z = ^n{z) has a single, simple solution z = within 
any such contour. Taking the infimum over the allowed e-values we obtain the estimate in 
f l4.34p . Moreover, using fl2.22p for g{z) = z we obtain the following explicit formula for the 
eigenvalues by an absolutely convergent series, 

°° 1 d™~^ 

= So{n) + e with e = a + E ^^^d^^^^)'" " 

m=l 

We now can write the equation in f l4.16p as 

«„(m) = g{0 , g{z) = Gn{z; m) (4.37) 

where g{z) fulfills the prerequisites of Lagrange's theorem 12.4.11 according to Corollary 14. 2. 41 
Thus f l2.22p gives an explicit formula for the an(m) by an absolutely convergent series. 
Moreover, inserting (14.340 in the estimate in (14.230 we obtain (I4.35p . 

To make this series more explicit we insert (12.230 and use the following multinomial 
series, 

r=0 eoA,---=0 llr=0 r- 

m 
1 

i=0 £o,---A=0 llr=0 r=0 

and similarly for g. This yields 

oo oo m— 1 r i£ 

^ = « + E E '5(Er=o'4,^)^(Er="M,m-l)(m-l)!n^ (4.39) 

m=l eoA,---/m-l=0 r=0 

and 

oo m— 1 oo 
m=l 1=0 £o,ii,...,ie=0 

x5(E'=i<, ^)(^ -mm- i)g^.e J] ^ • (4-40) 



E(^-a)' E ^(E^=o^.,^)5(E^=ir^r,^)-^nH'^ (4-38) 

1 ±r 



We observe that 



^r = ^^ma)-a)'^^^:\a). (4.41) 



r! da 

Recalling (I4.36P and inserting (14.150 we obtain (l4.3Up and (I4.3ip : we used 
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for complex parameters x^, which follows from the Leibniz rule. Similarly, recalling 

1 



9r 



r! da' 

and (14371) and using ( ]4.17p we obtain ( ]4.32p and ( 14.33p . 
The estimates in f l4.35p and Lemma I2.3.1I imply 



Gnia;m) = G^;\a;in) (4.43) 



|^n(x)| <5^|Mm)||/^(z)| <^(5(m,n) + gS.^-^(-^-"^)C') x 

m m 

C qi:,{K\m,\~Kjm,) ^ ^ ^- E, Kjn, (^Y., K\n,\ ^ ^ (^^^) 

for all b> 0, with C = 2B/[A - B - \a\ + ^{A-B - |a|)2 -4^2] and the constants in 
f l2.20p and fl4.20p . This proves that ^I/nlx) is uniformly bounded, and thus square integrable, 
on its domain [— 7r,7r]^. This completes the proof of (B). 

If we only know that a is such that fl4.18p holds true for some A > we cannot conclude 
anything about convergence. However, the results stated in (A) still hold true in the sense 
of formal power series; see Theorem 12.4.11 and Remark I2.4.2[ □ 



Remark 4.3.2 Since $n(^) is proportional to (see Remark 4-2.8 ) one can choose the 
parameter rj in lli4-36[ ) to be , and this shows that the m-th term in lli4-30[ ) and ^4-32^ is 
Qi^q^my Prom this we expect that these series have a finite radius of convergence in q^ in 
general. However, since the function (p{z) also depends on q^ , this is not easy to prove. 



Remark 4.3.3 The result in Theorem \4.3.1\ explains why the power series of in q 
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is very complicated: due to Ii3.11]) the term oc on the r.h.s in d^.i^p is 0{q'^^^~^); see 
Remark \C.2.1\ in the Appendix. Thus all such terms oc 7** with s = 2, . . . ,iN contribute 
to the power series coefficient oc of S^, and the complexity of this q^^-term therefore 
increases dramatically not only with i but also with N. 

Remark 4.3.4 Our proof of convergence relies on estimates in Section \4^ which are crude. 
In particular, we estimate the energy denominators [So{m) — So{n) — z] in l^.i^p and ^^4jJ^ 
by their smallest possible value which, if resonances exist, are assumed for resonances. 
However, resonances are rare, and it is easy to see that there always exists a constant 
Aq > such that \SQ{m) — So{n)\ > Aq for all m different from a resonance. It thus 
seemed to us that it should be easy to improve our estimates in Proposition \4.2.3\ enough 



to prove convergence of our series for arbitrary particle numbers N. However, despite of 
much effort (delaying the publication of this paper for two more years) we have not been 
able to do this up to now. Anyway, we have several independent reasons to believe in the 
existence of a finite radius of convergence for arbitrary parameters: (i) the results in IKTf . 
(a) the existence of variants of our solution method described in Section and Remark 2 



in Section\^ which both avoid resonances, (Hi) recent numerical results iBL^ . 
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4.4 A method to avoid resonances 



We now describe a generalization of our solution described in Sections 14. 1H4.3I which allows 
to avoid resonances. The idea is to use a variant of degenerate perturbation theory. 

We fix n G and consider the set of all corresponding resonances, i.e. all m G 
different from n such that —nj) = and So{m) = So{n). It is easy to see that there 

can be only a finite number of distinct resonances which we denote as n2, n^, . . ., n/j+i 
with R the number of resonances. We also set rii = n. For example, for = 2 and integer 
A we get -R = 1, 111 = (ni, and 112 = (722 — A, rii + A); see ( IS.lOp . 

We now apply Lemma 14.1.11 to (14. ip using the following projection, (PnQ;)(m) = a(m) 
for m = rij, J = 1, 2, . . . , i? + 1, and otherwise; A and B are as before. The general 
solution of (14. 6 p is now = J2j=i cjSnj with constants cj to be determined. This implies 
the following generalization of (I4.10p . 

n 

^„(x) = Y,Y. cjan,,(m)F„,(x) (4.44) 

J=l mez'v 



with an/(ni) obtained by setting n = nj in (14. lip and 



1 clef f if m = nj for J = 1, 2, . . . , i? 

1-1 



[[£^o(ni) - i^]]n [ i^oi™-) - £] otherwise 



(4.45) 



The equation determining the eigenvalues = ^o{'^)+^n and the coefficients cj is obtained 
from (14. 7p which can be written as follows, 

R+l 

SnCj = J2 '^r.j,nA^n)cK (4.46) 
K=l 

with 

*n.«W = -i;7' E tls^. . .rr ^. (4.47) 

s=i K=i Yll=i ^o{t^k + ELi ^e) - ^o(n) - z 

LL JJ n 

It should be possible to solve (I4.46P explicitly using a matrix version of Lagrange's reversion 
theorem. One now can expand about z = 0, and we expect that our estimates in Section 
can be made sharp enough to prove absolute convergence in a finite g-interval without 
restrictions on parameters. It would be interesting to explore this generalized solution 
algorithm in more detail, but this is beyond the scope of the present paper. 



5 Final remarks 

1. The functions J'n{z) defined in (ll.lOp can be expanded as follows, 

oo 

Mz) = J2JiA^)Q'' (5.1) 
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where i7o,n(z) are (essentially) the Jack polynomials |McDl [St] , as discussed in the intro- 
duction. It is interesting to note that for partitions n, all functions 

(^i^2---^iv)'Jn,Kz), £ = 0,1,2,... (5.2) 

are symmetric polynomials pTlj . It would be interesting to investigate if these polynomials 
have a combinatorial significance also for i > 1. 

2. Proposition 12.2.11 above suggest that the function in (12.11) has an expansion in eigen- 
functions ^'n(x) of the eCS Hamiltonian as follows, 

F(x;y) = 5^K„M/„(x)^;:M (5-3) 

n 

for some constants and the bar indicating complex conjugation, and our algorithm 
provides a means to extract from this the eigenfunctions. 

It is well-known that the eCS model has a family of independent (formally) self- 
adjoint differential operators of the form Hk = (— i)^ Sj=i ^+(1°^^^ order terms) for k = 
1,2, ... N which mutually commute, [Hk, Hi] = for all k,i = 1,2, ... , N, and including 
the total momentum operator P = Hi and the eCS Hamiltonian H = H2 in (11.11) |OPj . 
This suggests that the remarkable identity should be generalizable to all these differential 
operators. 

Conjecture: [Hk{x) - iffc(y)]F(x; y) = for all k (5.4) 

(for k = 1 the proof is trivial). Using this one could extend our results to the operator 
'^k=i ^kHk for arbitrary real coefficients hk. The resulting formulas would be similar to ours 
but with ^o(n) in fl2J3D replaced by the E5i=i {uj + j{N + 1 - 2j))\ The freedom to 
choose the parameters bj arbitrarily should allow to avoid the resonance problem completely, 
similarly as in |KTj . It would be interesting to find a direct proof of the identities in (15.41) 
and to derive and explore the explicit solution obtained with the operator Ylk=i ^kHk- 

3. To set our results in perspective we note that a formal series representation of eigen- 
functions and eigenvalues of the kind derived in this paper can be given for any quantum 
mechanical model. Indeed, assume that we want to diagonalize some self-adjoint Hilbert 
space operator H using some countable generating set /„ of the pertinent Hilbert space 
such that 

Hfn = EJn + (^«« = 0) (5-5) 

m 

for some constants En and Vnm- For example, if if is a Hamiltonian and a complete 
orthonormal basis then (/„, H fm) = En for m = n and Vnm otherwise, but what we 
describe here is equally true for any generating sets which are not orthonormal and/or are 
overcomplete. It is easy to see that our arguments in Section H] can be generalized to the 
present general case and imply that H has formal eigenfunctions 

i^n = 'Yan{m)fm (5.6) 
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with corresponding formal eigenvalues Sn = En + Sn where Sn is a solution of the equation 
Sn = with 




^ >' nK\ K1K2 ^ KsU 

K=i[[Ek - En- z]] 




where l/[[-Efc — z\\n = l/[Ek ~ z] for k ^ n and for /c = n. Moreover, the coefficients for 
these eigenfunctions are given by 



(this result is known as Brillouin-Wigner perturbation theory; see e.g. |Lo] and references 

therein). Using the formula in (12.221) one can obtain from this fully explicit formulas for 
the eigenvalues and eigenfunctions as formal series to all orders; see the formulas given in 
the proof of Theorem 14.1.21 This shows that it is not the existence of an explicit series 
solution which makes the eCS model special, but it is rather the details of this solution. To 
be more specific, the eCS model is special since there exists a basis /„ such that En and 
Vnm are given by simple explicit formulas, that Vnm only depends on n — m, and that the 
matrix Vnm is close to triangular which implies that the series solution converges (the latter 
property we only proved here in special cases, but we believe that it is true in general). 

It is worth mentioning that the formulas above shed some interesting light on quantum 
mechanical perturbation theory in general. In particular, they highlight that, for a given 
Hamiltonian H, it is the choice of the basis /„ which determines the usefulness of pertur- 
bation theory. Moreover, the formulas above hold true for any linear operator H and can 
be easily generalized to degenerate cases as described in Section 14.41 

4. Our algorithm was based on the remarkable identity in (12.51) which we first obtained 
using quantum field theory techniques \L1\ [L2] . While we latter found an elementary proof 
of this result (presented in Appendix \A.3\} . we feel that the quantum field theory proof is 
more illuminating since it not only shows that this result it true but also why. We also 
note that there exist interesting generalizations of this identity which we found using these 
quantum field theory results |L6j . 

5. To judge the usefulness of our results it is important study our series solution numerically. 
Numerical results for N = 2 were recently obtained by J.C Barba and show that the 
numerical convergence is much better than the estimates in Section 14.21 suggest [BL] . To be 
more specific, it seems that resonances are no problem in practice, and low order truncations 
of our series approximate the exact solution well not only for small g- values but up to values 
close to g = 1 and for a wide range of coupling parameters. 

6. For A = 1, g = 0, and partitions n the functions /n(z) in (12.91) are identical with the 
Schur polynomials |HLj . Our results therefore include simple explicit formulas for an elliptic 
deformation of the Schur polynomials which seem worth to be studied in more detail. 

7. We believe that the solution method presented in this paper can be extended to other 
examples, including the anharmonic oscillator and the Heun equation]^ We hope to come 

'''This was suggested to us by Vadim Kuznetsov. 




(5.8) 




25 



back to this in future work. 
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A Identities of elliptic functions 

In this Appendix we give an elementary proof of the Fact stated in fl2.1l) -( l?]2|) and on which 
our algorithm is based. For the convenience of the reader we also include the proofs of some 
properties of elliptic functions which we need. 



A.l Relation of V and p 

Here we prove (11. 2p and (11.61) . 

From the definition in (11.21) it is obvious that V{z), z G C, is doubly periodic with 
periods 2uJi = 27t and 2uj2 = i/5, it has a single pole of order 2 in each period-parallelogram, 
V{z) — is analytic in some neighborhood of 2; = and equal to cq in (II. 3p : we used 
q = exp(7ria;2/a;i) = exp(— /3/2). These facts imply the second identity in (11.21) (see e.g. 
|EMUT] . Section 13.12). 

To prove (II. 6p we note that 6{2z) equals, up to a constant, the Jacobi Theta function 
^i(^): 

^(^) = n i/4rr-\i ^ ^1(^/2) 

2r/^ nn=i(i-r") 

(see e.g. page 470 in [WW] ) . From the relation between -^i and the Weierstrass elliptic 
functions a, ( and p we therefore conclude (see e.g. page 473 in |WWj ) 

C(^) = -f log^(z) + ^ (A.2) 

and 

Piz) = -Ui^) = -f-2 log ^ W - - (A.3) 
dz dz"^ uoi 

where rji/uoi is a constant. To determine the latter constant we use the definition in (11.51) 
and compute 

00 

log e{z) = log[(z/2) - {z/2f/6] + log[(l - q'l' + g'" z'] + 0{z') = 

n=l 
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const. + log{z) 



' ^ oo 

24 ^ ^(T 

. n=l ^ 



J2n 



Recalling that pi^z) — z vanishes for z ^ one concludes from this and (1A.3P that 
Tji/uji = Co, recalling (11 .SI) . This together with (11.21) and (]A.3|) proves (11.61) . □ 



A.2 Proof of the identities in pTsn-pTl) 



We start with the following well-known identity for the Weierstrass elliptic functions ( and 

P, 

[Cix) + Ciy) + Ciz)]^ = Pix) + p{y) + p{z) ifx + y + z = (A.4) 
(this identity is given as an exercise on page 446 in [WW]). From flA.2p we conclude that 



(f){x) = 6'{x)/6{x) equals ({x) up to a term linear in x. Thus the identity in (IA.4p remains 
true if we replace ( by (j). This together with the trivial identity 



(j)ix)(f){y) + 0(a;)0(2) + 0(i/)0(^) = 
|[0(x) + <P{y) + <f>{z)]' - l[<f>{xr + <P{yr + <P{zy] 



implies (12.31) . 



□ 



A. 3 Proof of Lemma 12.1.11 



Let F = F(x; y) as defined in (12. ip . We compute 

d 



dx-i 



with 0(x) = 9'{x)/9{x), and thus 



^ = E ^^'^^^ -^k) A0(xj - yk) 



^ X'^cpixj - Xfc)0(xj - a;^) + ^ X'^(f){xj - yk)(l){xj - ye) 

-2 ^ \^(p{xj - Xk)(t){xj - ye) 



F . 



With that we compute straightforwardly 

F^Kdx] dyV 
which we write as a sum of four terms, W = Wi + W2 + W3 + W4, with 

j k^j 
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[x ^ y]' means the same terms but with the arguments Xj and Uj interchanged), 



= m 5Z ~ ^k)(t){xj - Xi) -[x ^y], 

and 

1^4 = ^ ^ - yk)(t){xj - yt) - 2A^0(xfc - xe)(j){xk - yj) - [x ^ y] = 

j,k e=/=k 

5Z 5Z ['^^'^(^i ~ yk)4>{xj - yi) + 2AV(yfc - yi)'P{yk - Xj) -[x ^y]. 

j,k e=/=k 

We first observe that the first two terms in W2 are invariant under x ^ y [note that 
(f)'{—x) = (f)'{x)], and therefore 

W2 = 0. 

We then write as follows [using 0(— x) = — 0(x)] 

W3 = ^ (-A^) (f){xk - Xj)(j){xj - Xi) + (j){xe - Xk)(t){xk - Xj) + 
j<k<e 



— \x ^ 



y] 



(j){xj - Xi)(j){xe - Xk) 
and using now the relation in fl2.3l) and f{—x) = f{x) we get 



j<k<i 



Finally, 



j k^j ij^j,k 

-{N - 2)x' J2 E [/(^^- - ^'^) - - 



j,k l=ik 



(piyi - yk)(t){xj - yi) -[x^y] 

where we wrote the same term in two different ways by renaming summation indices. We 
can now use the relation in fl2.3p again, and we obtain 

j,k e=/=k 
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where the terms even under [x y] canceled. Collecting all terms and using 
and 2/(x) = V{x) — + Cq we get 



+\^Nf{xj - Xk) -[x^y]-- 
5Z 5Z ['^'^'(^J ~ + - + 2X^f{xj - Xk) -[x ^ y] ■■ 

We thus see that WF is equal to the r.h.s. of fl2.2p . 



□ 



B Proof of Proposition 12.2.1 



We first observe that Lemma [2.1.11 remains true if we replace F(x;y) by 

F'(x; y) = ce'^^f=i(^^-^^) F(x; y) 



(B.l) 



for arbitrary constants P G M and c G C. [To see this, introduce center-of-mass co- 
ordinates X = Yl!j=i^j/^ x'j = {xj — xi) for j = 2,...,N, and similarly for the 
y's. Then H{x.) = —d'^/dX'^ + H^i'x'), and similarly for -ff(y). Invariance of (12. 2p under 
F exp [-iP(X - Y)N]F thus follows from + d/dY)F{x; y) = 0, and the latter 

is implied by the obvious invariance of -F(x; y) under xj xj + a, yj yj + a, a & M. 
Invariance of (12.21) under F — > cF is trivial, of course]. 

As shown below, one can find constants c and P such that the function -F'(x; y) above 
is the generating function for the -Fn(x) in (12. 8p as follows. 



F'{^-y) = j2h{^)^rUr^ ■■■i 



N ' 



6 = e'^^ 



nGZ 



with 



n+ = nj + ^{N +l-2j) 



and the series converges absolutely in the following region, 

1 < 161 < 161 < •■■ < l^jvl < q'^ 



(B.2) 



(B.3) 



(B.4) 



Moreover, in the same region one can change variables in the eCS Hamiltonian i/(y) from 
yj to 6 and expand 



[B.5) 



j<k vez 
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with the coefficients S^, in (12.141) . Inserting this in if(x)F'(x;y) = H [y) F' {'x; y) (imphed 

by Lemma (I2.1.ip ) and equating the coefficients of ^ '■■^at ^ 011 both sides we 

obtains the identity in (|2jId with £o(n) = Ej(^^)^ as in (|2A3|) . 

To prove (]B.2|) we note that the definition of the functions /n(z) in (12.90 is equivalent 
to the following generating function, 

Y[l<j<k<N '^i^i / ^k) _ \ - , x^_„j^_„2 



Moreover, using 

0{y) = |e'-/2r^/'e(0 for i = if |e| < 1 

following from (11.51) and the obvious identity sin(?//2) = |e"^/^e~'^/^(l — e'^), we find by a 
straightforward computation that -F(x; y) in ( 12. ip is equal to the expression on the l.h.s. of 
(IB. 6 II multiplied by ^E^ol^) in (II. 9p and the factor 

np— iA(j/j — ?/fe)/2 
e L-v- -n- - u~ i<fc = const e'^^T.'^=i^^j-y3)nQ-'^>^J:f^iiN+-^^'23)yjn 

This shows that (1B.2P holds true for certain constants P and c. We are left to prove 



V{y) = - ^ S^r'' for ^ = such that |^| < 1 (B.7) 



which obviously implies (lB.5p . For that we insert the identity 
4sin2[(y,±i£)/2] " 



PQ^^^^-^^ for all real and e > 



in ([L2D and obtain V{y) = - ^"^^^ v {i^"" + YZ=\ e"'"^(r + T'')] • Summing up the geo- 
metric series in m we obtain (iRTl) . 



C Further proofs 

C.l Proof of Lemma [27331 

We ffist prove the estimates in (12.180 using the definition in (12. lip : 

00 

|e(^)| < (1 + kl) n [(1 + ^'"^1^1) (1 + ^""/l^l)] = 

m=l 

and 

00 

|0(^)l > (1 - kl) n [(1 - ^''"l^l) (1 - ^""/I^D] = 

m,=l 
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for g2 < < 1. With that we can estimate the integrand in (12.91) as follows, 



"3 lij 



where the r.h.s. is constant on the integration paths in f l2.10p . Thus (12.91) implies the 
following upper bound, 



l/n(z)|<e 



n,>0(e- 



(C.l) 



for any choice of allowed parameters ej. To make the bound in (IC.ip more specific we 
choose 



P 



N + b' 



P 



1 + 26 



with b > arbitrary; note that these parameters satisfy the conditions in (I2.10p . We also 
observe that 

-Q(-x) < — < ' < Q(-q ) for < X < 1 

2^^ 1 + x l + g2 ^ ^ ^ ^ 



and 



> — — > — ^ — - > for < xo < X < Xi < 1 



1 — Xi 1— X I — Xq 

these estimates are easily proved using 0(±x)/(l =Fx) = n^i[l + 5'^" "F <f"'{x + 1/x)] and 
the fact that the function x + 1/x is monotonically decreasing in the interval < x < 1. 
Thus 



|/„(z)| < e^^"^^''^^-^!"^!) 



2\]N{N-l)X/2 



[(l-Xi)e(xo)/(l-Xo)r^ 

with Xo = miufc^nC"^* = e"^'-^"'"''-' and Xi = max^^ne"^'' = e"*^^^"^-*. Inserting e~^ = and 
pK = K we obtain Xq = q"^'"^^^ and Xi = q^~^ . □ 



C.2 Proof of the estimates in ( 14.251) 



In this appendix we derive the estimate (14.250 for the sum Ks{m) defined in (I4.24p . 
Representing the Kronecker delta by integrals, 



N 



and recalling (12.70 . y ■ = ^K^Vj^ — Vk^), we obtain 

N 



27r 



(C.2) 



a<j<k<N uel 
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for < El < 62 < ■ ■ ■ < Sn < P i^L the limit En I 0. The shifts in the integration paths, 
Vj ~^ Vj ~ ^^i' guarantee the absolute convergence of the sums in the integrand. 

It it important to note that, actually, this limit need not be taken: changing variables to 
= e'^J"'"'^^ the latter integral can be written as follows. 



with the integration paths Cj defined in fl2.10p . The integrand is analytic for all Ej obeying 
the condition in f l2.10p . as can be checked by recalling the definition of in (12.141) and 
= e~^. Thus Cauchy's theorem implies that integral on the r.h.s. in flC.2p is independent 
of the Ej. 

Inserting Ej = Je where < e < j3/N and using the triangle inequality and 5*1^ > we 
deduce from flC.2l) that 



j=i ~" / \j<k uei 

= Q-J:,j^rn, j |7| ^ ^ S.e-"'^^-^^' j . (C.3) 

\ j<k v&L ) 

To simplify this bound we use < |z/|gl'''~'^/(l — g^) following from the definition in fl2.14p . 
and thus 

1 ^-^(fc-J> _ 



j<k j<k u€Z ^ ^ 



^ ' j<k u=l 



JV-l oo 



^ ' 1=1 u=l 

1 / e"^^ Q-[f3-N6]-ee 



2\ 



[^^ ~ ^) (1 - e-fe)2 + \i - e-\P-m-'-^ f ) 



where we used = e ^ and changed summation variables. Inserting the estimate 

P 1 

< -7T if < p < 1 and < px^ < 1 



(1 — pxY (1 — x)2 
for p = exp{—[P — Ne]) and x = exp(— £:) we get 



^-1 -fe 



^ ^ (1 - g2) 2- (1 - e-^-)2 < 



N ^ e"^^ A^(A^-l)e- 
(W) tr (1 - ^"')' ^ (l-e-)3 
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since (1 — e '•^ ^^'^) < (1 — g^). To summarize, 



7|A^(A^ - l)e-^ 
(1 -e-=)3 




(C.4) 



Setting e = P/{N + b) for 6 > we get the bounds in (14.251) . 



□ 



Remark C.2.1 It is instructive to see in more details how the upper hounds for K ^{0) come 
about: It is easy to see that the contributions of lowest order in come from the identity 
in \3.11\) . for example for s < N the lowest order contributions are for u^. = — Ej^^^^-^ 
for r = 1,2,. . .,s - 1 and Us = Thus for s < N, Ks{0) = CN,s\l\"{S-iy~^Si + 

. . . = cn,s\'j\''^(i'^ + . . . with the dots indicating higher order terms, and, more generally, 
Ks{0) = cj\[^s\i\^q'^^^~^ + the constants cn^s are combinatorial factors. This and 

( [C^l j for e = (3/N and m = suggest the improved estimate in ^.29 ). 
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